Single hole and vortex excitations in the doped Rokhsar-Kivelson quantum dimer 

model on the triangular lattice 
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We consider the doped Rokhsar-Kivelson quantum dimer model on the triangular lattice with 
one mobile hole (monomer) at the Rokhsar-Kivelson point. The motion of the hole is described by 
two branches of excitations: the hole may either move with or without a trapped Z2 vortex (vison). 
We perform a study of the hole dispersion in the limit where the hole hopping amplitude is much 
smaller than the interdimer interaction. In this limit, the hole without vison moves freely and has 
a tight-binding spectrum. On the other hand, the hole with a trapped vison is strongly constrained 
due to interference effects and can only move via higher-order virtual processes. 

PACS numbers: 71.10.-w, 74.20.Mn, 05.50.+q 



I. INTRODUCTION 

Quantum dimer models in two dimensions are an ac- 
tive field of research due to their possible relevance 
to high-temperature superconductivity^ 2 ^ and to frus- 
trated magnetism,^ in the framework of the resonating- 
valence-bond (RVB) construction.— It was, however, re- 
alized that the RVB liquid state of the original Rokhsar- 
Kivelson (RK) dimer model on the square lattice is 
critical 6 and does not provide a stable phase. A search for 
a liquid with topological order 7 continued with quantum 
dimer models on nonbipartite lattices^ in particular on 
the triangular lattice. At a special value of the coupling 
constants (the RK point), the RK dimer model on the 
triangular lattice is proven to have exponentially decay- 
ing correlations jSilS and a gapped excitation spectrum is 
found numerically. 11 Furthermore, there is now convinc- 
ing numerical evidence that this gapped liquid phase is 
stable and extends within a finite parameter ranged 

In the RVB scenario of high-temperature superconduc- 
tivity, the superconductivity emerges upon doping a liq- 
uid of singlet bonds with charge carriers i 5 i 13 i 14 i 27 Very 
early it was therefore proposed to study quan tum dimer 
models containing monomers^ (see Refs.[TlllTlfl9||20| 
for more recent studies). When analyzing a hole in the 
background of dimer or spin liquids, one should take 
into account vortexlike excitations inherent in topologi- 
cal liquidsi2J£ These topological excitations were dubbed 
visons in the context of the corresponding Z2 gauge 
theory. 21 There is now evidence that visons constitute 
the gapped excitations of the undoped triangular-lattice 
quantum dimer liquid . n > 12 It is known that visons can 
bind to holes and change their statistics^*^ 

In this work, we study the properties of a single-hole 
excitation in the case of a doped triangular-lattice quan- 
tum dimer model at the RK point and illustrate the 
hole-vison binding. We find two branches of excitations: 
one for the hole itself and the other for a hole-vison 
bound state. The energy-momentum dispersion for both 
branches of excitations is computed in the perturbative 
regime of small hole hopping. The combined effects of lat- 



tice frustration and Z2 flux lead to a quadruply degener- 
ate vison-hole branch and reduced bandwidth. These re- 
sults may have interesting implications for RVB physics. 

This paper is organized in the following way. In Sec.HTl 
we introduce the model and show the existence of two 
branches of excitations. In Sec. Mil we calculate per- 
turbatively the dispersion of the non- vison branch of the 
hole excitation. In Sec. HVi the dispersion of the vison 
branch is calculated. Finally, in Sec. El we summarize 
and discuss our results. 



II. DOPED DIMER MODEL, TOPOLOGICAL 
SECTORS, AND TWO BRANCHES OF 
EXCITATIONS 



We consider the quantum dimer model on the trian- 
gular lattice doped with mobile holes. We choose the 
simplest form of the hole-hopping term which involves 
rearrangement of one dimer. The Hamiltonian reads 
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RK+hole — 

E H (l//>cri + h - c -) + v (I ~)cri + I//X//I)] 
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+ E[- s d")(A|+h.c.) + 2 M |/.)(/.|] 

A 

= H RK + H s , (1) 

where the first sum is performed over all three orienta- 
tions of rhombi and the second sum is over both up and 
down triangles and over all three possible positions of the 
hole on the triangle. 

We consider the model [Eq. ((TJ] at the RK point, t = 
V = 1, in the sector with a single hole. At s — u > 0, the 
Hamiltonian has the usual "supersymmetric" properties 
of the RK point: its ground state is exactly known and 
given by the equal-amplitude superposition of all possible 
states) 2 " and the quantum mechanics can be mapped onto 
a classical stochastic dynamics in imaginary time.— We 
further consider the hole term H s as a perturbation in 
s <C 1, u <C 1. To simplify the formulae, we assume 
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u = s > 0, but our results are trivially extendable to 
u ^ s. 

In the unperturbed Hamiltonian Hrk , the position of 
the hole x is a static parameter. We consider the hole 
on the infinite lattice (or, equivalently, on a large finite 
lattice far from the boundary) . In such a setup there are 
two degenerate ground states of Hrk for each hole posi- 
tion. They correspond to two disconnected (topological) 
sectors H ± (x) of the Hilbert space, characterized by the 
values ±1 of the vison operator 

\r( \ / i \No. of dimcrs intersecting T x /o\ 



for some contour T x connecting the hole position x to 
infinity (in a finite system to the boundary) .— ^ The cor- 
responding ground states are given by the sums over 
all dimer coverings in the respective topological sector 
and are denoted as ip^(x). Note that while the label- 
ing ± of the two sectors H ± (x) depends on the choice of 
the contour T x , the sectors themselves do not. Chang- 
ing the contour T x amounts to possible interchanges 
H + (x) <-> TL~{x) and, therefore, iPq(x) «-> ^Pq(x). This 
ambiguity reflects the Z 2 degree of freedom in labeling the 
topological sectors, and it will play an important role in 
the motion of the hole with a trapped vison. Technically, 
this Z2 gauge may be fixed by specifying (arbitrarily), 
for each 2, a reference dimer covering which belongs to 
H+(x). 

The two topological sectors ~H ± (x) differ by the par- 
ity of the dimer intersection at infinity and hence are 
indistinguishable by any local operator (since all corre- 
lation functions are short-ranged in the RK model on 
the triangular lattice^^). Therefore, for excitations ob- 
tained from the ground states by local operators (in 
the vicinity of x), one can establish a one-to-one lin- 
ear mapping between the states in H + (x) and in TL~{x). 
Taking odd and even combinations of the corresponding 
states ip + (x)±tp~ (x), we obtain the decomposition of the 
Hilbert space into even and odd sectors H e '°(x). Those 
even and odd sectors correspond to the non-vison and 
vison sectors of excitations, respectively, introduced in 
Ref. [111. 

The key observation for our discussion is that the 
Hamiltonian [Eq. {T])] preserves the decomposition into 
H e '°(x) at every point x. While it is obviously true for 
Hrk and the potential part of H Sl one can also easily 
check that the hopping part of H s does not have matrix 
elements between H e (x) and H.°(x') for neighboring sites 
x and x' . Hence, the excitations of the moving hole can 
also be classified into two branches: the non-vison branch 
[contained in ® x Ti e {x)] and the vison branch [contained 
in ^^"(i)]. This splitting into even (non-vison) and 
odd (vison) branches is a generic feature of any pertur- 
bative mixing of topological sectors in quantum dimer 
models. 




FIG. 1: These two types of hopping processes set opposite 
correspondences between the sectors Ti ± (x) and Tl ± (x'), and 
therefore cancel each other in Eq. (|6]). 

III. NON-VISON BRANCH 

The energy spectrum of the non-vison branch can be 
easily calculated to first order in the perturbative expan- 
sion in H s . We can fix the phases of all RK ground states 
ipg(x) = ipQ (x) + ipQ (x) by taking the linear combination 
of all dimer coverings with the amplitude one (up to nor- 
malization). Then, the problem of the moving hole maps 
onto the tight-binding model with the hopping amplitude 

h = ~(%(x)\H s \%(x')) (3) 

for nearest-neighbor x and x' (here and in the following 
we always assume normalized states). This amplitude 
may be converted into an expectation value in the RK 
model with a static hole, 

t 1 =2s{Mx)\~)(~\Mx))=2sw> ( 4 ) 

iVi 

where N\ and N3 are the numbers of dimer coverings with 
one site and one three-site triangle removed, respectively. 
The ratio N3/N1 is well defined in the limit of the infinite 
system and can be computed numerically with a suitable 
method. We have calculated this coefficient with a Monte 
Carlo simulation similar to that in Refs.|T(3 and[ll| (using 
clusters of toroidal geometry with up to 17x17 sites), 
with the result iV 3 /./Vi = 0.229 ± 0.001. 

Taking into account the potential term in H s and per- 
forming the Fourier transformation in x, the dispersion 
of the hole without a vison takes the form 

E\s_ = —2ti (cos fci + cos k 2 + cos fc 3 — 3) , (5) 

where k\, k 2 , and fc 3 are the projections of the vector k 
on the three lattice directions (with k\ + k 2 + k% = 0) . 

IV. VISON BRANCH 

The hopping of a hole with a trapped vison is more 
complicated. The phases of the odd-sector ground states, 
V'o( :c ) = i^oix) ~ *Po( x )> cannot be synchronized invari- 
antly for all x, which reflects the frustration of the vison 
motion~i The freedom of the Z 2 gauge [the contours T x 
in Eq. ^ or, equivalently, the reference dimer configura- 
tion for each x] corresponds to the choice of the overall 
sign for the states in Tt°(x). 
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FIG. 2: The two trajectories of the hole exactly cancel each 
other in the second order of the perturbation theory [Eq. (0] 
for (a) nearest-neighbor and (b) next-nearest-neighbor hop- 
pings. The nonvanishing second-order terms are (c) on-site 
and (d) next-next-nearest-neighbor hoppings. 



Regardless of the chosen Z 2 gauge, the hopping ampli- 
tude vanishes to first order, 

(ro(x)\H s \r (x')) = o (6) 

for nearest-neighbors x and x' . This can be seen as the 
cancellation of the two types of hopping processes from x 
to x 1 , corresponding to two possible dimer flips (Fig. 1). 
Each of those dimer flips maps each of Ti.^ (x) into one of 
W ± (a;'). The change in topological sector depends on the 
chosen gauge, but the correspondence between the two 
sectors 7i (x) and the two sectors W ± (cc') is opposite for 
the two types of flips^ As a result, the corresponding- 
processes connecting two ground states ipo(x) and ipo(x') 
exactly cancel each other. 

A nontrivial hopping appears only to higher order in 
perturbation theory for some trajectories. The second- 
order hopping amplitude 

*a= E ^mx)\H s \rjx')){r n {x')\H s \r {x')) 

(7) 

involves excitations ip^(x') of Hrk with energies E n . 

Similarly to the cancellation of the nearest-neighbor 
hopping amplitude to first order in perturbation, one 
can show the cancellation to second order of the hop- 
ping processes x — > x' — > x" connecting nearest-neighbor 
and next-nearest-neighbor sites [processes (a) and (b) in 
Fig. [2]. One can verify that, in those cases, processes 
symmetric with respect to the line xx" exactly cancel 
each other. 

The only nontrivial hopping in the second order occurs 
for trajectories x" involving two links along the 

same direction [i.e., for the on-site energy correction and 
for the next-next-nearest-neighbor hopping, processes (c) 
and (d) in Fig. |5J. The corresponding next-next-nearest- 
neighbor hopping amplitude [Fig. [^d)] to second order 




- 2 4 2 4 2 




FIG. 3: (a) We fix the relative gauge at the next-next- 
nearest-neighbor sites by relating the reference configurations 
in TL + (x) and TL + (x") via two consecutive dimer flips on oppo- 
site sides of the line xx'x" . (b) The four sublattices connected 
by the hopping of the vison-hole bound state. This composite 
excitation can only hop by multiples of two lattice periods. 

in perturbation [Eq. |(7J)] may be expressed via dynamic 
correlation functions in the RK model with a static hole 
at position x' , 

/>OC 

h = / dt{%(x)\H s e- H * Kt H s \%(x")) 
Jo 

/>OG 

= s 2 / dtl(t), (8) 
Jo 

where 

I(t) = (Mx')\Px X >e~ HRKt Px>x"\Mx')) (9) 

and 

P*x> = \/.)(/.\-\\')(\'\ 

p X ' X » = I V)( VI - 1.\ XA I- (io) 

The dynamic correlation function I(t) is well defined in 
the limit of infinite system size and does not depend on 
the topological sector in this limit. It may be computed 
with a classical Monte Carlo method as in Ref.hjl Us- 
ing clusters of toroidal geometry and up to 17x17 sites, 
we find f^dt I{t) = -1.51 ± 0.08 (observe that it is 
negative). 

Note that the sign of t% in Eq. ||5J)] corresponds to a 
particular relative gauge choice at points x and x": the 
reference dimer coverings at x and x" are connected by 
two dimer flips on opposite sides of the line xx'x" (Fig. [3]). 
One can show that this local gauge convention for any 
two sites separated by two lattice periods can be con- 
sistently extended to a global gauge on the sublattice of 
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such sites (with the period of this sublattice equals twice 
that of the original lattice). There are four such sublat- 
tices (Fig. |3|), and the hole-vison excitation hops on each 
of them independently, without a possibility to cross over 
to another sublattice. The resulting dispersion relation is 
that of the tight-binding model with the doubled lattice 
constant and the hopping amplitude given by Eq. ([B"|. 

= -2t 2 (cos2fci + cos2fc 2 +cos2fc 3 ) + s . (11) 

The on-site energy £o is equal to that in the non-vison 
sector in Eq. ([5]). To leading order in s, it is given by 

Eq = 6t X . 

The hole-vison excitations with dispersion (TTTj) are 
quadruply degenerate (by sublattice) for each value of 
k in the Brillouin zone of the doubled lattice. While we 
have explicitly demonstrated this degeneracy to second 
order, it can be extended to all orders of perturbation 
theory. In fact, this degeneracy is determined by the 
symmetries of the original Hamiltonian [Eq. ((T|)] in the vi- 
son sector and can thus be promoted from a perturbative 
argument to the exact spectrum. The exact degeneracy 
can be proven using the translational invariance of the 
Hamiltonian, together with the symmetry under point 
inversion (rotation by tt) and time reversal. Physically, 
this degeneracy can be understood as the cancellation of 
virtual processes for the flux-carrying excitation on the 
frustrated triangular lattice. 

Finally, let us note that, while our derivation of the 
vison-hole spectrum was formally done at the RK-point, 
its form and degeneracy are the same in the whole liquid 
phase away from the RK point (estimated to extend to 
the region 0.8 < j < 1 in Ref. [H), provided the hole 
hopping is small. Only the numerical coefficients in the 
hopping amplitudes t\ and <2 get modified in this case. 
Furthermore, our results equally apply when more than 
one hole is present in the system as long as the holes are 
sufficiently far apart and do not interact with each other. 

V. SUMMARY 

In this paper, we have calculated the dispersion of a 
single mobile hole in the RVB liquid phase of the doped 



RK quantum dimer model on the triangular lattice. We 
find two branches of excitations: one for the bare hole 
and the other for a hole-vison bound state. The effective 
motion of the hole-vison state is strongly modified by the 
Z 2 flux associated with the vison. Interference effects due 
to lattice frustration reduce the bandwidth of this type 
of excitation and lead to additional degeneracies. These 
are general properties, which should be observed in any 
doped Z2 RVB liquid on frustrated lattices. 

In our specific model [Eq. {TJ], the energy of a static 
(s = 0) vison-hole bound state equals that of a hole with- 
out a vison. In other words, the vison does not cost any 
energy if placed in a hole (while in the bulk, its energy is 
a finite fraction of t, see Ref. [TT1 ). In the limit of a small 
hopping amplitude s, the energy of a static excitation is 
split, with the bandwidth proportional to s for the bare 
hole and to s 2 /t for the hole-vison bound state. As a 
result, the two branches intersect each other, with the 
minimum of energy (the ground state) corresponding to 
the hole without a vison. For some k in a region close to 
the boundary of the Brillouin zone, the vison-hole bound 
state is lower in energy than the bare hole. In a more 
general quantum dimer model (or in other RVB-type sys- 
tems), however, one may imagine the situation where the 
hole-vison bound state constitutes the ground state (in 
our dimer model, this may be achieved, for example, by 
adding ring exchange of dimers around a hole). In such 
a case, the doped holes spontaneously generate visons, 
which, in turn, may lead to further interesting effects, 
e.g., the modification of the statistics of holes.— ^ 
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